Composite fermions in a long-range random magnetic field: 
Quantum Hall effect versus Shubnikov-de Haas oscillations 



A. D. Mirlin, 1 '* D. G. Polyakov, 2 <t and P. Wolfle 1 

1 Institut fur Theorie der Kondensierten Materie, Universitat Karlsruhe, 76128 Karlsruhe, Germany 
2 Institut fur Theoretische Physik, Universitat zu Koln, 50937 Koln, Germany 

We study transport in a smooth random magnetic field, with emphasis on composite fermions 
(CF) near half-filling of the Landau level. When either the amplitude of the magnetic field fluctu- 
ations or its mean value B is large enough, the transport is of percolating nature. While at B = 
the percolation effects enhance the conductivity a xx , increasing B (which corresponds to moving 
away from half-filling for the CF problem) leads to a sharp falloff of a xx and, consequently, to the 
quantum localization of CFs. We demonstrate that the localization is a crucial factor in the inter- 
play between the Shubnikov-de Haas and quantum Hall oscillations, and point out that the latter 
are dominant in the CF metal. 



PACS numbers: 73.40.Hm, 71.10.Pm 

The idea of composite fermions (CF) has been re- 
markably successful in accounting for the transport prop- 
erties of strongly correlated two-dimensional electrons in 
high magnetic fields, especially near half-filling v = \ of 
the lowest Landau level. Within this model, each elec- 
tron is replaced by a fermion carrying two flux quanta of a 
fictitious magnetic field oriented oppositely to the exter- 
nal field B ext . At half- filling, the average of this Chern- 
Simons field cancels the external one exactly, thus offer- 
ing an elegant explanation of the observed gapless state 
with a well-defined Fermi surface. Away from v — A, the 
composite fermions experience, in the mean-field sense, 
a net effective field B = (1 — 2^)B ex t. An integer num- 
ber p of filled Landau levels of the fermions then corre- 
sponds to the fractional quantum Hall (QH) plateau at 
v = pj (2p ± 1). This mapping of the singularities in the 
strongly interacting electron system at fractional fillings 
onto integer QH transitions is a particularly attractive 
feature of the CF picture. 

Confidence that the CF picture does indeed capture 
the essential physics of the v — ^ problem is strength- 
ened by a number of observations of Fermi-surface fea- 
tures at half-filling in experiments with acoustic waves 
and in direct transport measurements jqj. In partic- 
ular, in the vicinity of half-filling the dissipative resistiv- 
ity p xx (B ext ) exhibits magnetooscillations || which look 
very much like the familiar Shubnikov-de Haas (SdH) ef- 
fect. However, there is an essential difference: the os- 
cillations near v = \ appear to be much more strongly 
damped as compared to the conventional metallic phase. 
Generally speaking, this is in line with the fact that p xx 
at B = is orders of magnitude larger than at B ext — 0. 
More specifically, the enhanced damping is understood 
in terms of strong scattering on spatial fluctuations of 
the Chern-Simons field. The fluctuations of the effective 
magnetic field B(r) = [1 — 2v(r)]B ext around the mean 
value B are generated by the inhomogeneous distribution 
of the electron density and therefore the local filling fac- 
tor v(r) due to screening of the impurity potential. This 



mechanism of scattering dominates [|J the conductivity 
of the CF metal, whereas the direct scattering on the 
scalar random potential is of little importance. 

This Letter is concerned with the disorder-induced lo- 
calization of the composite fermions and its implications 
for transport measurements. In particular, we address 
the question as to the nature of the magnetooscillations 
near half-filling. Though it has become customary to 
identify the oscillations with the SdH effect H, we ar- 
gue that they are due to the quantum localization re- 
sponsible for the QH effect. Specifically, the SdH and 
QH oscillations coexist with each other; what we show is 
that the SdH oscillations are much more strongly damped 
as compared to those induced by the quantum localiza- 
tion. We also demonstrate that the effective magnetic 
field creates a classical percolation network, which dras- 
tically suppresses the CF conductivity. It is this classical 
localization effect that enhances the quantum localiza- 
tion leading to the QH oscillations. 

The model we deal with is that of non-interacting 
fermions in a random magnetic field (RMF) B + B(r) 
with the average B and the correlator (B(0)B(r)) — 
B^F(r). We wish to calculate the CF conductivity as 
a function of B. The measurable electron resistivity pf, v 
can be directly expressed in terms of the CF one [0: 
Pxx = Pxx, p% y — Pxy — 2h/e 2 . Disorder is assumed to be 
created by ionized impurities randomly distributed with 
a sheet density rii (for simplicity let it be equal to the 
electron density n) in a thin layer separated from the 
electron gas by an undoped spacer of width d rT 1 !" 1 . 
The correlation function of the RMF is then parameter- 
ized by B = k F /V2d and F{r) = (1 + r 2 /4d 2 )~ 3 / 2 , 
where k F = Ann and K = c = e = 1. 

The RMF is characterized by two length scales: the 
correlation radius d and the cyclotron radius in the field 
Bo, Ro = hp /Bo. Defining the parameter a = d/Ro, we 
can distinguish the weak- RMF regime a <C 1, where the 
mean free path I ^> Rq ^> d, and the regime of strong 
fluctuations a>l, where one should expect drastic de- 
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viations from the Drude picture. Though a = l/\2 for 
the CFs at half-filling, it is instructive to get a feeling for 
the problem by exploring these two limiting cases. Also, 
in principle one can change a by varying the ratio of the 
impurity and electron densities, since a = (rii/2n) 1 / 2 . 

We start with the limit of weak B where the qua- 
siclassical treatment is accurate Q. At a <C 1, CF 
trajectories are only slightly bent on the scale of d, so 
that the Born approximation |^^| is valid. Accord- 
ingly, for the transport scattering time one gets l/r tr = 
(Bl/mkp) J™drF(r) = 2a 2 v F /d, where the CF effec- 
tive mass m = Uf/vf is introduced. The conductivity at 
zero B then reads (in units of e 2 /h) 

a xx = kpd/Aa 2 , a«l. (1) 

If the size of inhomogeneities had been strictly zero 
("white noise" disorder), one would have obtained at the 
quasiclassical level the magnetoresistance Ap xx /p xx — 0. 
To get a finite magnetoresistance, one has to include 
the effect of the field B on the "ballistic" scale d. To 
this end we write the equation of motion in the form 
d(j>/dt = oj c + B[r(t)]/m, where cj) is the angle of the ve- 
locity on the Fermi surface and oj c = B/m. The velocity 
correlation function vp (cos\<j)(t) — 0(0)]) is then given by 

v F (cos (— J *dtiS[r(fi)]jy Here (...) denotes an en- 
semble average over exact classical trajectories r(t). In 
the limit of small a, the classical dynamics of the system 
is fully chaotic (uncorrelated diffusion) and so Z?[r(t)] can 
be treated as a Gaussian variable with respect to the dis- 
order averaging. Accordingly, the conductivity is 

a xx +ia xy =^- rdte^SM 
2m J 

S V = o^a [ dt 4 dt2 (B[r(h)]B[r(t 2 )]) . (2) 
zm Jo Jo 

The double integral can be rewritten to give S(t) = 
(B /m) 2 J*dt'(t - t') (F(\r(t') - r(0)|)). The argument 
of F is readily expanded at t <C r tr in powers of to c : 

|r(t)-r(0)| ~ « F i[I - M) 2 /24] . (3) 

We can identify two different contributions to the mag- 
netoresistance. One comes from a B dependent correc- 
tion to r^ 1 = lim^oo t~ 1 S(t). Substituting Eq. (3) into 
S(t) we get Tt r l (B) ~ T t ~ 1 (0)[l + (d 2 /2R 2 c ) ln(Z_/d)], where 
R c = vf/uJc is the Larmor radius in the field B. This cor- 
rection would yield a positive Ap xx / p xx . Next, we should 
take into account that the relaxation kernel is not exactly 
the simple exponential e~ l l Ttr \ specifically, Sit) cx t 2 at 
t < d/vp. This gives a 4 times larger negative contribu- 
tion to the magnetoresistance, so that in total 

A £xx= _3_f_ I 
p xx 2Rj. d 



This result agrees with that derived in Ref. 7 by a differ- 
ent method. The sign of Ap xx /p xx found for a <C 1 is 
opposite to the one observed experimentally near v = i. 
The discrepancy is likely due to the fact that the condi- 
tion of weak RMF, I 3> d, is not met in the experiments. 

Let us now turn to the strong- RMF regime, a ^> 1, 
keeping B = 0. The seemingly innocent assumption 
about the chaotic character of the particle dynamics, 
which enabled us to represent the conductivity in the 
form (2), is not valid any more. Most particles are now 
out of play since they are caught in cyclotron orbits drift- 
ing along the closed lines of constant B(r) ("van Alfven 
drift"). In the adiabatic limit, their drift trajectories 
are periodic and so do not contribute to the conductiv- 
ity. Still, however large Bq is, there are classical paths 
which are not localized and percolate through the sys- 
tem by meandering around the lines of zero B(r). The 
conductivity is determined by the particles that move 
along these extended "snake states" || . Note that there 
is one single percolating path on the manifold of the 
B(r) — contours; yet, the conductivity is nonzero since 
the snake-state trajectories form a bundle of finite width, 
R s ~ d/a 1 / 2 . The conducting network is made up of 
those snake-states that can crossover from one critical 
zeio-B line of length L s ~ a 7 ^ 6 d (for a review of the per- 
colation theory see Ref. 9) to another. The coupling of 
two adjacent percolating clusters occurs near the critical 
saddle-points of B(r), which are nodes of the transport 
network. The characteristic distance between the nodes, 
i.e. the size of the elementary cell £ s , is then a 2 ^ 3 d B. On 
length scales longer than £ s , the particle dynamics can be 
viewed as fully stochastic. We estimate the macroscopic 
diffusion coefficient as D ~ v s F) 8 , where v s ~ L S R S /^ 2 is 
the fraction of particles residing in the delocalized snake- 
states, D s ~ £ 2 vf/L s their diffusion coefficient. We thus 
have D ~ vfR s and, correspondingly, 

o xx — k F d/a 1/2 , a > 1 . (5) 

It is worth noting that by comparison with the Born ap- 
proximation [Eq. (1)], the conductivity is ~ o?l 2 times 
larger (though the localization effects are strong and 
naively one might have expected the opposite). This is 
consistent with the experimental observation jq] that the 
CF conductivity at half-filling is a factor of ~ 5 larger 
than the perturbative-in-a result Q a xx = kpd/2. Let 
us also note that a xx given by Eq. 5 is larger by a factor 
of ~ a 1 / 2 than that obtained for a ^> 1 in Q by using an 
"eikonal approach" . The fault in M is not with the qua- 
siclassical approximation itself, but with the method of 
disorder averaging, which is incorrect in principle at large 
a; in particular, it neglects the localization of particles 
and the percolating character of the transport through 
the snake states. Whether the experimentally observed 
positive magnetoresistance can be understood within the 
percolation picture is not clear at present. 
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We now consider the regime of strong B, which is de- 
fined for a ~ 1 by the condition B 3> -Bo? or equiva- 
lently, R c <C d. Since d/R c — kpd/2p, where p is the 
number of filled CF Landau levels, and experimentally 
kpd ~ 10 ~ 15, the range of p where the magnetooscilla- 
tions are observed, p < 7, undoubtedly requires a strong- 
B treatment. At B ^> So, the CF dynamics is a slow 
van-Alfven drift of the cyclotron orbits along the lines 
of constant B(r), i.e. the conductivity is again deter- 
mined by a percolation network of trajectories close to 
the B(r) = lines. Naively one may well think the per- 
colation picture is very much like the one we dealt with 
at zero B and a 3> 1. In actual fact, there is a crucial 
difference. Specifically, now there is no stochastic mix- 
ing at the nodes of the percolation network: unlike the 
snake states at B = 0, the rapidly rotating cyclotron or- 
bits pass harmlessly through the critical saddle-points of 
B(r) without changing to the adjacent cell. In the high- 
B limit, the mixing occurs on the links of the network 
and is only due to the weak scattering between the drift 
trajectories. 

In order to calculate the conductivity at B ^> Bq, we 
first need to integrate out the fast cyclotron rotation, tak- 
ing care not to lose the effect of the non-adiabatic mixing. 
A similar question for the electron system in a random 
scalar potential was recently addressed in ]l0| . We write 
the equation of motion in the form 

z(t) = iv F fdt'e^o #"W>] ) (6) 
Jo 

where z(t) — x + iy stands for the exact trajectory r(t) = 
(x,y) in the field B + B(z), and — uj c + B(z)/m. 

It is convenient to introduce the guiding center coordi- 
nate ((t) = z(t) - i? c [C(i)]e lv(t) , where R c is the local cy- 
clotron radius at the point of the drift trajectory and 
ip(t) = J" d£'fi c [£(£')] is expressed in terms of the local 
cyclotron frequency. The velocity of the guiding center is 

given by C(i) =iv F U$a d *' n M*')]_ e M*)) _ j?^. After 
averaging over the cyclotron motion (denoted by (. . .) c ), 
the first term leads to the usual expression for the drift 



velocity v d = (C) c 



mv 2 F /2(B + B) 3 \ [VB x (B + B; 

At this level the drift occurs strictly along the lines of con- 
stant B(r), so that the conductivity would be zero. The 
time dependence of the cyclotron radius in the second 
term gives rise to the leading non-adiabatic contribution 
Svd = ((2vj/vF)e' lip )c, allowing the particle to perform 
a random walk around the B(r) = contour. Because 
of the rapidly oscillating factor e lip , Sv d will be exponen- 
tially small. We define the diffusion coefficient associated 
with the motion perpendicular to the drift trajectory as 
D x = Ilim^ 0o t- 1 ([/ Sv d (t)dt] 2 ) d cx e- w , where (...)«, 
denotes the average over the impurity ensemble along the 
drift lines. 

Expressing the phase ip by integration along the con- 
tour, tp — lo c j dlvj (I), and shifting the integration into 



the complex plane one observes that the exponent W is 
determined by the phase tp picked up at the singular point 
of the correlator (B(Q)B(r)) at r = 2id. Consequently, 



W 



In 



i Uc f* td dlv-i(l) 



(7) 



where the integration should be done along the straight 
line connecting the points I — and I — 2id. It follows 
that the dominant contribution to D± comes from rare 
fluctuations in which the drift velocity greatly exceeds 
the typical value vf(Bq/ B) 2 . Indeed, W can be written 
as a sum of two terms, 



Wl = - 



i r d 2 k 



</.r k 



(2tt) 2 (v dx v dx ) 



W 2 



IU> C 



lid 



dx 



(8) 



Here W\ represents the probability for the optimum fluc- 
tuation u2 (r) to occur, Wi stands for the probability of 
the non-adiabatic scattering on this fluctuation. Solu- 
tion of the variational equation SW/8v^ x — then yields 
v dx( r ) = v F(Bo/B)J : (r/d), where J 7 is a dimensionless 
function. We thus get 



W = c{B/B^ 2 



(9) 



where the numerical coefficient c ~ 1. Choosing v dx {r) oc 
(v dx {0)v d x(r)) as a trial function with the variational pa- 
rameter v^ x (0), we obtain the estimate c ~ 1.6. 

The non-adiabatic mixing of drift trajectories yields 
the conductivity <j xx ~ kpd, where 5 is an effective width 
of the links of the percolation network. This 6 obeys the 
equation S 2 ~ D±L(5)/v d , which is the condition of con- 
nectivity of the network. Here L(S) ~ d(d/S) 7 ^ 3 is the 
characteristic perimeter of the cells || . Consequently, 



Vxx = k F d x /(— ) exp 



3c / B x 1 



13 V B, 



(10) 



where f(x) is a power-law function. The conductivity 
is seen to fall off sharply beyond the scale B ^ Bq. 
This unusual manifestation of classical localization is also 
found in non-interacting systems with a scalar random 
potential Jl(]]. Finally, it is worth noting that, however 
sophisticated the dissipative conductivity network may 
be, the macroscopic Hall conductivity a xy in the metal- 
lic system at B > Bo assumes the collisionless form 
<j xy — 2-rm/B = p > a xx JTTJ . 

1 1 a v ; r «, found the monotonic component of a xx (B), we 
turn to the magnetooscillations. These are convention- 
ally interpreted as the SdH effect for the CFs. We want 
to compute the disorder-induced Dingle factor of the 
SdH oscillations Sa^ H cx cos(4n 2 n/B) exp[-S r (B)]. At 
B >• B the damping exponent is given by S r = | ( ( I >2 ), 
where $ is the random- B flux through the cyclotron or- 
bit in the field B 0. Had the CF dynamics been fully 
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chaotic, (. . .) here would have meant a disorder aver- 
age over all possible configurations of B(r), which yields 
S r = 4ir 2 p 4 /(kFd) 2 The result might look puzzling, 
since this S r is numerically by far larger Q than what is 
observed experimentally However, in the limit B 3> Bq, 
only the particles that drift along the percolating tra- 
jectories with zero B{y) contribute to a xx . Accordingly 
whereas the above S r describes the damping of the total 
density of states, the averaging we need to do should be 
performed along the contours B(r) — 0. The result is 



S r = (2l7r 2 /4)p 8 /(/fc F d) e 



fill 



Both the extremely sharp falloff of the amplitude of the 
oscillations with increasing p and the numerical value of 
S r are in fair agreement with the experimental data avail- 
able to date ||; yet, there are good reasons to question 
this picture. The point is that the CF system, like a 
conventional metal, will exhibit magnetooscillations even 
if those of the density of states are totally neglected. 
Specifically, in the diffusive regime, the quantum inter- 
ference of scattered waves leads to the QH oscillations 
g§ Sag H cx cos(2ira xy )eiq)[-Si(B)], where 



Si = 2lT(7 xx 



(12) 



The confusion that one might have at this point is over 
the notion that the QH effect transforms into the SdH os- 
cillations with decreasing magnetic field, which is man- 
ifestly not true. Likewise, let us stress that the obser- 
vation of the QH oscillations at p > 1 does not re- 
quire exponentially low temperatures, - while the char- 
acteristic localization length does grow exponentially fast 
with p, the oscillations come from the interference on the 
small scale of the effective mean free path. Note that 
the periods of the QH and SdH oscillations coincide at 
B ^S> Bq, but the damping factors are very different. The 
quasiclassical approach which led us to the exponential 
falloff of a xx with decreasing p fails and the QH effect 
shows up at pQ H ~ k F d/ ln 1 ' 2 (/ci?(i), where a xx drops to 
unity. At p < p® H the oscillations are due to the QH 
effect. This picture is in good agreement with experi- 
ment: from the typical p xx ~ 0.02 -j- 0.03 at p ~ 6 -f- 7, 
where the first resistivity minimum is observed, one finds 
a xx ~ p xx p 2 — 1 -j- 1.5. Furthermore, at smaller p < 5, 
the values of a xx at the maxima are close to i, thus con- 
firming that the CFs are indeed in the QH regime. As 
for the SdH oscillations, Eq. (11) tells us that they might 
only be observable at p < p SdH ~ (k F d) 3 ^ < p QH . The 
last inequality means that in the limit k F d 3> 1 there is 
no room for the SdH effect. In practical terms, however, 
pQH anc j pSdH are numer i ca Uy close at the experimen- 
tally relevant k F d ~ 10 4- 15, so that the SdH effect still 
can contribute to the first couple of minima (p = 7,6). 

In conclusion, we have studied the fermion kinetics 
in a smoothly varying RMF with mean B. We calcu- 



lated the conductivity in the regime of strong RMF and 
zero B, when the transport is determined by percolating 
snake-states. We demonstrated that increasing B leads 
to the classical localization of fermions, the key signa- 
ture of which is the exponentially sharp drop in a xx . We 
showed that this yields a strong enhancement of the mag- 
netooscillations, thus explaining why they are observed 
at large p, where the SdH oscillations are negligible. We 
analyzed the interplay between the SdH and QH oscil- 
lations and argued that the latter are dominant in the 
limit kpd ^> 1. 
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